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The energy distributions of the Gamow- Teller strength are studied for even-even Xe isotopes with 
mass numbers from 124 to 142. A self-consistent microscopic formalism is used to generate the single 
particle basis, using a deformed Skyrme Hartree-Fock mean field with pairing correlations in BCS 
approximation. The Gamow- Teller (GT) transitions are obtained within a quasiparticle random 
• phase approximation (QRPA) approach using a residual spin-isospin interaction in the particle-hole 

and particle-particle channels. We then discuss the pairing BCS treatment and the determination 
of the ph and pp residual interaction coupling constants. We study the GT + and GT~ strength 
£Nj ' distributions for the equilibrium nuclear shapes, which are an essential information for studies of 

charge-exchange reactions and double-/3 processes involving these isotopes. 
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. I. INTRODUCTION 

\Q , The isotope 136 Xe has been recently used as a moving target colliding with a hydrogen gas jet in a first test to 
proof the feasibility of the EXL experimental approach at FAIR-GSI [1]. Charge-exchange (p,n) reactions on this Xe 
isotope were consequently measured. Although this kind of facilities are intended to explore highly unstable nuclei, 
some stable isotopes like those under study here are normally used as initial test targets. Therefore, having reliable 
information on their nuclear structure turns out to be of primary importance. From 136 Xe, we extend our study to 
other stable isotopes, 124_134 Xe, and also to the neutron-rich region 138 > 140 > 142 Xe. The latter undergo /3~-decay and 
| their half-lives provide us with another piece of experimental information. 

Gamow- Teller transition matrix elements can be extracted from the measured forward-angle charge- exchange data 
[2]. At high incident energies and at forward angles, the nuclear states are probed at small momentum transfer. 
Therefore, only the central parts of the isovector effective interaction contribute to the cross-section. Furthermore, 
because of the small momentum transfer, a multipole expansion leads to a simple relation between the measured 0° 
cross sections and the corresponding allowed /3-decay transition rates (L=0). Assuming that even taking into account 
the projectile distortion effects in DWBA the transition amplitudes would still approximately factorizc into a nuclear 
reaction part and a nuclear structure part, the 0° charge-exchange cross sections are proportional to the corresponding 



Fermi and Gamow- Teller matrix elements. The proportionality factor is parametrized in terms of a distortion factor, 
a kinematic factor and a volume integral of the effective interaction [2] . 

Theoretical GT ± strength distributions, as the ones discussed here, can therefore be used to predict cross sections of 
various charge- exchange reactions under the appropriate kinematic conditions. Some of the charge-exchange reactions 
corresponding to GT~ processes are: 

136 Xe(p,n) 136 Cs ; 136 Xe( 3 He, 3 H) 136 Cs ; 136 Xe(d,2n) 136 Cs ; 

while some examples of those corresponding to GT + transitions are: 

136 Xe(n,p) 136 I ; 136 Xe( 3 H, 3 He) 136 I ; 136 Xe(d,2p) 136 I. 



In addition, the Xe isotopes are of special interest regarding double (3 processes [3] . Different nuclear models have 
been developed by many groups to calculate quantitatively the double (3 matrix elements, as described in recent review 
articles [4]. Table I shows all the transitions of this type involving Xe isotopes as parents or daughters. The GT^ 
strength distributions which will be obtained here can be used to calculate the transition amplitudes of the initial and 
final ground states going to the virtual QRPA-excited states of the intermediate odd-odd nucleus. After computing 
the overlap between the intermediate states coming from the parent and from the daughter nuclei, the two-neutrino 
double beta decay matrix elements as well as the half-lives can be calculated. Such a calculation has indeed been 
carried out in [8] for 128 Xe, 130 Xe and 136 Xe, which are j3~ / j3~ parents or daughters. 
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Finally, Xe isotopes are of considerable interest because they belong to a typical shape transitional region [9], in 
which there are experimental indications of triaxial deformation in some isotopes [10]. Present theoretical triaxial 
calculations are mainly based on algebraic models [11]. Nevertheless, in this work we assume axial symmetry. As we 
shall see, for the deformed Xe isotopes with two equilibrium shapes there are no critical changes in the Gamow- Teller 
strengths at these two shapes, and a similar behavior may be expected when considering possible triaxial shapes. In 
any case, axial deformation is a crucial ingredient of the formalism that gives rise to new features in the Gamow- Teller 
strength distributions different from those obtained within a spherical treatment. 

Although some of the isotopes studied may be spherical, in this work we are not restricted to those but we deal with 
a large number of Xe isotopes whose equilibrium shapes are unknown. For most of the isotopes considered here, there 
is no clear experimental evidence of whether they are spherical or deformed. Under these circumstances, a deformed 
approach is always preferable over a spherical one because a deformed formalism contains the spherical shape as a 
particular solution. 

We have also found in the past [12,13] that the GT strength distributions may depend, in some cases significantly, 
on the deformation of the decaying nucleus. But we notice that this dependence has to be studied case by case since 
it is very sensitive to the fragmentation and crossing of levels generated by the deformation. We think it is indeed 
worth studying the degree of sensitivity of the GT strength of these Xe isotopes to deformation. 

Following this introduction we present in Sec. II a brief description of the theoretical framework. Sec. Ill includes 
our results regarding HF+BCS energies and GT^ strength distributions, together with a discussion on the pairing 
treatment. Finally, Sec. IV contains the main conclusions of our work. 

II. THEORETICAL FRAMEWORK 

We describe here briefly the theoretical formalism used, whose details can be found in ref. [12]. We carry out a 
deformed Hartree-Fock calculation with the effective nucleon-nucleon density-dependent Skyrme interaction Sk3 [14], 
assuming axial deformation and time reversal symmetry. The single-particle wave functions are expanded in terms of 
the eigenstates of an axially symmetric harmonic oscillator in cylindrical coordinates using eleven major shells. Pairing 
correlations between like nucleons are included in BCS approximation either taking fixed pairing gap parameters (A w 
for protons and A„ for neutrons) or taking fixed pairing strength parameters (G ff and G v respectively). We refer to 
these two types of calculations as HF(Sk3)+BCS(A) and HF(Sk3)+BCS(G), respectively. They yield single particle 
energies and wave functions together with their occupation probabilities for protons and neutrons separately. A 
quadrupole constrained HF+BCS calculation [15], where the intrinsic quadrupole moment is constrained, is also 
performed to obtain the deformation dependence of the ground state energy. 

The pairing energy gaps A are determined phcnomenologically, and the pairing strengths G are obtained from them 
in an indirect way, as will be described later on. Within BCS approximation, both parameters are related by the so 
called gap equation: 

A = Gj2wi, (2-1) 

i 

where Vi y u, are occupation and non-occupation probability amplitudes of the i th single particle level subject to the 
condition vf + uf = 1. In order to determine the value of G that reproduces a given value of A, one should take into 
account that this depends on the active energy range and number of levels considered. In our case, we include all HF 
single particle sates in our basis above and below Fermi level. It is important to stress that the occupation probability 
amplitudes are computed in each iteration of the HF method, and are used to calculate the one-body density and 
mean field of the next iteration, so that one gets new single-particle wave functions, energies and occupation numbers 
at each iteration. Therefore, the selfconsistcnt determination of the binding energy and deformation includes pairing 
correlations from the beginning. 

To describe Gamow- Teller excitations we add to the quasiparticle mean field a separable spin-isospin residual 
interaction in the particle-hole (ph) and particle-particle (pp) channels, which is treated within QRPA. The advantage 
of using separable forces is that the QRPA energy eigenvalue problem is reduced to find the roots of an algebraic 
equation. The ph part 

V& = 2*gr E (-1)^+/^, (3+=J2(u\a K \n)ata^ (2.2) 

AT=0,±1 TTV 

is responsible for the position and structure of the GT resonance [12,16-18]. The corresponding coupling constant, 
Xph 1S obtained in a consistent way from the same energy density functional as the Hartree-Fock mean field through a 
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second derivative with respect to the nucleonic density and averaging the contact interaction over the nuclear volume 
[12]. The pp part consists of a proton-neutron pairing force, which we introduce as a separable force [17,19,20], 
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(2.3) 



The coupling constant n^J may in principle be derived consistently with the HFB or HF+BCS mean field through 
a second derivative with respect to the pairing tensor of the energy density functional. This derivation would be 
analogous to the way in which the ph force is obtained as the second derivative with respect to the density. In our 
theoretical scheme the proton-neutron pairing interaction is neglected in the construction of the mean field to avoid 
mixing of even-even and odd-odd isotopes in the intrinsic state. Only pairing between like particles is included in the 
BCS approximation. This implies that the particle-particle interaction in the proton-neutron channel is in principle 
undetermined. Therefore the coupling constant k^J is fitted to the phenomenology, as for example to reproduce 
half- lives as it is usually done [17,18]. 

The pnQRPA phonon operator for GT excitations in even-even nuclei is written as 



(2.4) 



where ir and v stand for proton and neutron, respectively, a + (a) are quasiparticle creation (annihilation) operators, 
ujk are the RPA excitation energies, and X%{f , the forward and backward amplitudes, respectively. It satisfies 



|0) = 0; 



?t K |0> = \u>k) ■ 



(2.5) 



when acting on the QRPA ground state of the parent nucleus |0). 

The technical details to solve the QRPA equations are well described in Rcfs. [12,16,19]. Here we only mention 
that, because of the use of separable residual forces, the solutions of the QRPA equations are found by solving first a 
dispersion relation, which is an algebraic equation of fourth order in the excitation energy ui. Then, for each value of 
the energy, the GT transition amplitudes in the intrinsic frame connecting the ground state |0) to one phonon states 
in the daughter nucleus \ujk), are determined by using the normalization conditions of the phonon amplitudes. They 
are given by 



where 
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with 



q-KV UjyV^Zj^ J q^ V — VjyU-^Zj^ 

It is a simple matter to find out that the Ikeda sum rule 



£ [(Af-^) 2 - {M+ K ) 2 J = 3(A - Z) 



(2.8) 



(2.9) 



holds in RPA approximation provided all the eigenvalues contained in the basis space are included in the sum, so that 
the orthonormalization conditions arc satisfied. In practice, the strength functions are calculated up to an energy 
u> < Scut) where E cut is such that Ikeda sum rule is fulfilled up to a few per thousand discrepancy. Typical energies 
used in our calculations are E cnt = 30 MeV. The number of configurations involved in this mass region for this energy 
range is typically over one thousand. 

Once the intrinsic amplitudes in Eq. (2.6) are calculated, the Gamow- Teller strength B{GT) in the laboratory 
frame for a transition IiKi(0 + 0) — ► IfKf(l + K) can be obtained as 



B ± {GT) = K J / M / 1^1 7 * M *)| 2 = { S *tfi {4>k s \Pt\4>of + 26 Kf ,i {4>k s \pf\ 4>of} ■ 

Mi,M f ,n 



(2.10) 



To obtain this expression we have used the initial and final states in the laboratory frame expressed in terms of the 
intrinsic states \4>k >, using the Bohr and Mottelson factorization [21]. Theoretical (3~ -decay half- lives are calculated 
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by summing up all the energetically allowed transition probabilities in (2.10), in units of g^/Air, weighted with phase 
space factors, up to states in the daughter nucleus with excitation energies below the Qp- window. 

One may wonder whether, in deformed nuclei, the calculated GT strengths may contain spurious contributions from 
higher angular momentum components in the initial and final wave functions. As mentioned above, the GT strengths 
are calculated in the laboratory frame in the factorization approximation of Bohr and Mottelson. Using angular 
momentum projection techniques [22], the angular momentum projection can be carried out through an expansion 
in inverse powers of the angular momentum operator component perpendicular to the symmetry axis < J\ >. This 
expansion, to lower order, provides a factorization approximation formally identical to that of Bohr and Mottelson. 
Thus, the effect of angular momentum projection is to a large extent taken into account. An upper bound to 
contributions from higher angular momentum components is proportional to < J\ >~ 2 , with values of < J\ > 
ranging from 10 to 40 in the case of the deformed Xe isotopes. Therefore, exact angular momentum projection in 
deformed Xe isotopes would lead in all cases to less than a few percent effect in the GT strengths. In the cases where 
the shape is spherical, there are not high angular momentum contributions to the GT strengths. 

It may be questioned whether it is correct to introduce additional BCS parameters and residual interactions on top of 
the Skyrme interaction, which is already an effective interaction. Indeed, if the effective force used in constructing the 
mean field were to be the most general possible interaction, one should not include additional parameters. However, 
this is not the case for the Skyrme interaction. The parameters of the Skyrme force are determined by requiring 
that they reproduce the nuclear compressibility, as well as the total binding energies and charge radii of magic 
nuclei in spherical selfconsistent calculations. It is well known that the effective Skyrme interaction and its existent 
parametrizations are suitable to generate the optimal HF mean field of spherical and deformed nuclei [14,23]. The 
particle and hole eigenstates of the mean field, which form the canonical basis, have highly nontrivial wave functions 
that contain a mixture of many harmonic oscillator shells, when expanded into a harmonic oscillator basis. In 
addition to deformation, for non-closed shell nuclei, one has to take into account the pairing correlation effects which 
are important when the level density around the Fermi level becomes large. 

It is also well known that the effective Skyrme interaction is not suitable to generate the quasiparticle mean field, 
since in the fit of the Skyrme forces no attention is paid to the realistic character of the pairing matrix elements, and 
extensions of the Skyrme HF method have been developed over the years. To include pairing correlations in the mean 
field one possibility is to do Hartree-Fock-Bogoliubov calculations using cither finite range forces, like the Gogny force, 
or contact density dependent pairing interactions [24]. All of them are extensions of the Skyrme forces specifically 
designed for this purpose. The other possibility is to perform BCS calculations in the canonical basis using either 
phcnomcnological fixed gap parameters or fixed pairing strengths, as originally proposed by Vauthcrin. This is the 
path followed in our paper. This path has been proved to be successful to study the properties of ground state and 
low spin excited states in open shell nuclei [12,14,23,25]. 

Concerning the residual interactions, we mentioned above that the same effective Skyrme interaction is used to 
generate a separable particle-hole residual interaction. The separable interaction simplifies enormously the calculation 
and still contains the main characteristics of the contact force. The quasiparticle energy density functional obtained 
with the effective Skyrme and pairing interactions that we used does not contain any dependence on particle-particle 
interactions in the proton-neutron channel. Therefore we have to introduce a proton-neutron particle-particle residual 
interaction in the usual way as a separable force with a coupling strength that we fit to the measured half lives. We 
would like to recall that a bridge between Skyrme HF and RPA calculations for excited states was established long 
time ago [26] by using a particle-hole force in the RPA, which is determined by the second derivatives of the HF energy 
with respect to the density. But there is no guarantee that the derived force is good for excited states. This is so 
because Skyrme forces are constructed for the description of ground state properties. Finally, we would like to notice 
that similar schemes based on Skyrme HF+BCS+RPA have been frequently used in the literature with successful 
results [27]. 

III. RESULTS AND DISCUSSION 

A. Equilibrium deformations 

A deformed Hartrec-Fock mean field calculation is performed using a Sk3 Skyrme force with a constraint of the 
quadrupole deformation given by the parameter (3 = Q p /(Z{r 2 )) 1 where Q p is the proton quadrupole moment 

and (r 2 ) is the charge mean square radius. A pairing interaction between like nucleons within BCS approximation is 
also included, keeping fixed the pairing energy gap (A) or the pairing strength (G). From this calculation we obtain 
the ground state energy as a function of the quadrupole deformation (3. 
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In principle, when a pairing force with fixed strength G^.j, is used, the pairing energy gaps A n>u depend on 
the strength of the pairing interaction as well as on the occupation amplitudes of the single-particle states (see 
Eq.(2.1)). This last condition amounts to say that pairing energy gaps are deformation-dependent. Actually, we 
obtain these gaps first phenomenologically from the odd-even mass differences by means of a symmetric five term 
formula involving experimental binding energies [7], and keep them fixed to carry out a deformation-constrained 
HF+BCS(A) calculation. Next we fix the pairing strength so as to reproduce the gap parameters at the deformation 
of the ground state, and we carry out a deformation-constrained HF+BCS(G) calculation. In this way we perform 
a fixed pairing strength calculation, which is conceptually more appealing, but still profiting from the experimental 
information available for the gaps. 

Fig. 1 shows the HF+BCS energies for 124 ~ 142 Xe, using Sk3 Skyrme interaction. The pairing interaction is included 
in both fixed gap (dashed line) and fixed strength (solid line) treatments, each of them in a different curve whose 
absolute minima have been separated 1 MeV for a better comparison. For the same isotope, both curves show energy 
minima at very similar deformations, the ones corresponding to the fixed pairing strength treatment being slightly 
smaller in absolute value. One can also observe that the energy barrier at the spherical region of 124 ~ 132 Xe and of 
140-142-^-g j g j esg p ronounC ed when the pairing strength is fixed, and in this case both minima have a very similar 
energy (the prolate one being generally the ground state, except for 126 Xe). As the mass number increases from 
A=124 to A=138, the deformations of the equilibrium shapes decrease and eventually converge to a spherical shape. 
The two final isotopes show only one equilibrium deformation in the prolate region. 

Underneath each of these energy-deformation graphs we plot the corresponding pairing gaps as a function of the 
quadrupolc deformation for the fixed pairing strength calculation. The vertical lines join the ground states from fixed 
gap treatment with the corresponding pairing gaps at this deformation from the fixed strength treatment. These 
values are those coming from the aforementioned odd-even experimental mass differences. 

In our calculations with fixed G values the binding energies show a tendency to increase at those deformations 
where the pairing gaps reach a minimum. This looks contradictory since one may expect the opposite because the 
smaller the pairing gap is, the lower is the contribution of the pairing energy to the total binding energy. However, 
it indicates that the minima of the pairing gaps appear at similar deformation to that where the volume and the 
spin-orbit term contributions to the binding energy are maximum. On the other hand, comparing calculations with 
G fixed to those with A fixed, one sees in 1 that for (3 values where A takes larger values the binding energy increases 
more compared to that obtained with lower (A fixed) value. 

In Table II we show the quadrupole deformation (3 of the equilibrium shapes according to the HF(Sk3)+BCS(A) and 
HF(Sk3)+BCS(G) calculations. We also show for comparison the results from independent theoretical calculations 
obtained from selfconsistent relativistic calculations [29] as well as from phenomenological nonrelativistic calculations 
[30]. These results also indicate the existence of deformed solutions, which agree with those obtained here. Upper 
limits of the ground state deformation obtained from experimental B(E, 2^ — > 0^) transitions are also included [28]. 
Table III shows the pairing gaps, obtained from experimental binding energies, and the pairing strengths reproducing 
these gaps at the ground state deformation. 

B. Gamow- Teller strength distributions 

The spin-isospin residual interactions in the particle-hole and particle-particle channels are treated here within a 
quasiparticle random phase approximation. The particle-hole residual interaction coupling constant Xph is obtained 
consistently with the Hartrec-Fock mean field, and their values vary from 0.21 MeV in 124 Xe to 0.19 MeV in 142 Xe. 
We have used an average value for all the isotopes under study, \ph = 0.2 MeV. In the case of the particle-particle 
residual interaction, we have chosen the value of the coupling constant k pp ho as to reproduce the experimental half- 
lives of the three unstable Xe isotopes: 138 Xe (T 1/2 =844.8 s), 140 Xe (T 1/2 =13.60 s) and 142 Xe (T 1/2 =1.22 s). A good 
agreement between calculated and experimental (3~ half-lives for the three isotopes is reached with k pp =0.07 MeV, 
provided we use the standard attenuation factor 0.77 for spin matrix elements as in previous works. 

The single particle energies and occupation probabilities at the equilibrium nuclear shapes are obtained from a 
HF(Sk3)+BCS(G) calculation. Fig. 2 shows GT~ strengths (in g 2 A /(A-K) units) as a function of the excitation energy 
of the daughter nucleus after the transition. Discrete and gaussian-folded distributions are shown, the latter being 
more suited to compare with experimental results regarding the GT strengths themselves or the cross sections of 
charge-exchange reactions obtained from them. The range of the excitation energy from to 30 MeV includes the 
resonance, which appears at around 13 MeV in 124 Xe and moves slightly toward higher excitation energies as the 
mass number increases, reaching 25 MeV in 142 Xe. 

Fig. 3 shows the same calculations but for GT + transitions. As expected from the Ikeda sum rule, Eq. (2.9), the 
scale of the strengths is much smaller in this case. Table IV contains the summed GT ± strengths, their difference, 
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the value of 3 (-/V — Z) and the fulfillment of the Ikeda sum rule in percentage for the prolate shape of every isotope 
(which is generally the ground state), and for the oblate shape when there is a second minimum. The Ikeda sum rule 
is fulfilled up to a very high degree of accuracy in all the cases. 

The fragmentation observed in the GT + strength is reduced as the prolate energy minimum moves to the spherical 
region ( 134 . 136 438j^ e ^ ^ double peak structure becomes then apparent, which was responsible for the larger width of 
these resonances in comparison with the ones in GT~ distributions. A second peak with less strength, in the near 
spherical shapes, appears 8 MeV below the biggest one and moves accordingly with it as the mass number changes. 
It is worth noticing the large single peak which appears at very low excitation energies in 126 Xe and also in 124 Xe, 
reaching in the prolate shape of this last case a strength of 0.34 g 2 A j(Ai:\ as indicated in the figure. This strength 
corresponds to a dominant GT transition from a K 7 ' = 9/2+ proton state to a K 7 ' = 7/2+ neutron state connecting 
the proton g 9 / 2 shell with the neutron g 7 / 2 shell. The occupation probability of the neutron state is small enough to 
allow the transition in 124 Xe and 126 Xe, but when the number of neutrons increases, this state becomes blocked for 
GT transitions. 

From the GT 1 * 1 profiles of 124_132 Xe it is obvious that no clear distinction can be made between oblate and prolate 
deformations. It is only possible to distinguish both deformations in some cases when small energy windows are 
explored, as for example in the low B(GT + ) energy window of 124,126 Xe. Similar studies on the effect of deformation 
in the GT strength distributions were done in the neutron-deficient Hg-Pb-Po region [31], and in the A~70 mass 
region [12]. 

The effect of deformation on the GT strength distributions can be observed more clearly in Fig. 4, where we compare 
spherical and deformed QRPA results. We show the examples of 128 Xe, where two equilibrium shapes, oblate and 
prolate, are obtained, and 140 Xe, where a prolate shape is predicted. In Fig. 4 plots upward correspond to deformed 
calculations with prolate shapes, while plots downward correspond to spherical calculations. The left panels show the 
GT~ strengths and the right ones the GT + strengths. As we can see in Fig. 4 the main effect of deformation is the 
stronger fragmentation of the strength, which is particularly clear on the GT + strength distributions because of the 
smaller scale. The positions of the peaks are also changed from spherical to deformed in a different way for each case. 

In order to compare the GT strength distributions obtained here with those coming from a HF(Sk3)+BCS(A) 
calculation of the single particle levels, we show in Fig. 5 the distributions corresponding to both BCS pairing 
treatments (fixed gap and fixed strength) for the prolate and oblate equilibrium shapes of 128 Xe. The results are very 
similar, as could be expected given the fact that the equilibrium nuclear shapes occur at very similar deformations 
and binding energies in both BCS treatments (see Fig. 1). In particular, for the prolate shape both distributions are 
almost identical as expected because the values of the G n>l , parameters were chosen to reproduce the A„. j!y parameters 
precisely at this deformation. It is interesting to know that also at the oblate minimum these pairing strengths nearly 
reproduce the values of the phenomenological pairing gaps, since as seen in Fig. 1 they are very close to the ones 
at the prolate minimum. Therefore the distributions for the oblate case with both pairing treatments are also very 
similar. In the case of spherical equilibrium shapes, as for example 136 Xe, the GT strength distributions from fixed 
pairing gap and fixed pairing strength calculations are indistinguishable. 

IV. SUMMARY AND CONCLUSIONS 

The Xe isotopes are of considerable theoretical interest because they participate in a variety of double-beta decay 
processes and because they belong to a nuclear shape transitional region. In addition, from the experimental point 
of view, the stable Xe isotopes have been used as moving targets in charge-exchange reactions to test new facilities, 
where the unstable Xe isotopes will be explored in the near future. The present work has addressed these topics by 
predicting stable nuclear shapes and GT strength distributions, which are a fundamental tool to calculate single and 
double beta transition matrix elements and half-lives, as well as cross sections of charge-exchange reactions under the 
appropriate kinematic conditions. 

For even-even Xe isotopes with mass numbers from 124 to 142, we have studied the GT strength distributions 
using a deformed pnQRPA formalism with ph and pp spin-isospin separable residual interactions. The quasiparticle 
mean field is obtained from an axially deformed HF approach, with the Skyrme interaction Sk3, including pairing 
correlations between like-nucleons in BCS approximation using either fixed gaps (A^) or fixed pairing interaction 
strength {G VjV ). The HF+BCS mean field has been also used to consistently determine the ph coupling constant for 
every isotope, whose average value has been finally used for all of them. The pp coupling constant has been fixed to 
approximately reproduce the half-life of the three /^"-unstable Xe isotopes included in this work. 

From the energy-deformation curves, an oblate-prolate shape coexistence is predicted in 124 ~ 132 Xe with a low 
energy barrier between them. 134 ~ 138 Xe are predicted to be spherical, whereas the prolate shape seems to be strongly 
favored in 140_142 Xe. In general, a fixed pairing strength calculation increases the binding energy of the spherical shape 
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region. The deformation dependence of the pairing energy gaps from a fixed pairing strength calculation has been 
also shown. The self-consistent quadrupolc deformations of the ground state derived within the HF+BCS procedure 
are in agreement with independent theoretical calculations as well as with the experimental upper bounds extracted 
from B(E2) transitions in the whole chain of Xe isotopes under study. 

The GT - strength distributions for the equilibrium shapes are dominated by a single peak moving to higher 
excitation energies and gathering more strength as the number of neutrons increases, as expected. In the cases where 
two equilibrium shapes are predicted, there is no strong dependence of the GT strength distribution on the equilibrium 
shape, at least when a wide range of excitation energy is considered. 

In the case of the GT + transitions, the strength is more fragmented, giving rise to a richer structure in the energy 
distribution. This fragmentation decreases as the energy minima move to the spherical region, where a double-peaked 
resonance appears as observed in 134 ~ 138 Xe. The lightest isotopes, 124 Xe and 126 Xe, show a very high peak from a 
single transition at very low excitation energies, but this transition is blocked in the isotopes with higher number of 
neutrons. With such a complex structure, the influence of the sign of the nuclear deformation (oblate or prolate) in 
124_132 Xe is more apparent on the GT + strength distributions than on the GT - transitions, but in any case it does 
not seem to be critical. However, the GT strength distributions obtained from spherical or deformed shapes show 
different features related to the energy location of the main peaks and to the fragmentation of the strength. From the 
accumulated GT 1 * 1 strengths up to 30 MeV of excitation energy, it has been shown that the Ikeda sum rule is fulfilled 
up to a very high percentage (see table IV). 

Theoretical studies of GT strengths and related observables as, in particular, charge-exchange reaction cross sections 
are necessary to help in the event simulation work of these kind of processes, that will be measured at the new FAIR- 
GSI facility. Theoretical work on this direction is in progress. 
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TABLE I. Double beta processes involving Xe isotopes, with their experimental half-lives [5,6] and Q-values from experi- 
mental masses [7]. When a /3 + //3 + transition is indicated, f3 + /EC and EC/EC are also allowed (with Q-values 1.022 MeV 



and 2.044 MeV lower than the 


one shown, respectively). 






Transition 


2/?-process 


T1/2 (yr) exp. 


Q- value [MeV] 


Xe — > Te 


(3 + /(3 + 


> 2.0 x 10 [5] 


2.866 


_a_o — > ie 






U.oy / 


128 Te ^ 128 Xe 


/T//T 


= 2.5 ± 0.3 x 10 24 [6] 


0.867 


130 Te ^ 130 Xe 


(3~/p- 


= 0.9 ±0.1 x 10 21 [6] 


2.529 


130 Ba - 130 Xe 




> 4.0 x 10 21 [5] 


2.610 


132 Ba - 132 Xe 


EC/EC 


> 3.0 x 10 20 [5] 


0.840 


134 Xe - 134 Ba 




> 1.1 x 10 16 [5] 


0.830 


136 Xe - 136 Ba 




8.1 x 10 20 [5] 


2.468 



TABLE II. Quadrupole deformation (3 of the 124_142 Xe equilibrium shapes from a HF(Sk3)+BCS calculation obtained with 
fixed pairing gaps A as well as with fixed pairing strengths G. Results from Refs. [29] and [30] are also given for comparison. 
Also given are the experimental values obtained from B(E2) transitions [28]. 



A 


Pth. 


(A fixed) 


Ptn. 


(G fixed) 


13th. [29] 


(3th. [30] 


\(3exp. | 




prolate 


oblate 


prolate 


oblate 






from B(E, 2+ -» 0+) [28] 


124 


0.24 


-0.19 


0.22 


-0.17 


0.215 


0.208 


0.264 (8) 


126 


0.19 


-0.18 


0.18 


-0.15 


0.186 


0.170 


0.1881 (30) 


128 


0.16 


-0.16 


0.15 


-0.12 


0.160 


0.143 


0.1837 (49) 


130 


0.13 


-0.13 


0.11 


-0.10 


0.128 


-0.113 


0.169 (6) 


132 


0.11 


-0.10 


0.07 


-0.10 


-0.070 


0.000 


0.1409 (46) 


134 


0.05 


-0.05 


0.01 




0.000 


0.000 


0.120 (10) 


136 


0.00 




0.00 




-0.001 


0.000 


0.086 (19) 


138 


0.03 




0.01 




-0.002 


0.000 


0.0309 (18) 


140 


0.15 


-0.10 


0.12 




0.104 


0.116 


0.1136 (25) 


142 


0.17 


-0.12 


0.16 




0.141 


0.145 





TABLE III. Pairing parameters for the HF+BCS calculation in 124 142 Xe. Pairing gaps A„ j7r are obtained from experimental 



binding energies [7]. The pairing strengths G v ,tt are those reproducing these gaps at the ground state deformation, which is 
also indicated (see text for details). 



A 


A, [MeV] 


A w [MeV] 


(3th. gs 


G v [MeV] 


G* [MeV] 


124 


1.32 


1.35 


0.24 


0.114 


0.132 


126 


1.31 


1.33 


0.19 


0.111 


0.133 


128 


1.27 


1.32 


0.16 


0.108 


0.130 


130 


1.25 


1.31 


0.13 


0.108 


0.129 


132 


1.18 


1.24 


0.11 


0.107 


0.125 


134 


1.01 


1.12 


0.05 


0.107 


0.119 


136 


1.44 


0.98 


0.00 


0.121 


0.112 


138 


1.00 


1.20 


0.03 


0.098 


0.120 


140 


0.96 


1.06 


0.15 


0.086 


0.144 


142 


1.03 


1.06 


0.17 


0.084 


0.116 
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TABLE IV. Calculated summed GT~ and GT + strengths (in <?a/(47t) units) for the ground state and the first + excited 
state (in brackets) of 124 ~ 142 Xe, from a HF(Sk3)+BCS(G) calculation. The difference between both summed strengths is 
compared with the value of 3(JV — Z) to check the fulfillment of the Ikeda sum rule (in percentage). 



A 




EGT+ 


T.GT" - EGT+ 


3(iV - Z) 


% 


124 


48.62 


1.15 


47.47 


48 


98.90 




(48.34) 


(0.84) 


(47.50) 




(98.96) 


126 


54.33 


0.78 


53.55 


54 


99.17 




(54.23) 


(0.71) 


(53.52) 




(99.11) 


128 


60.16 


0.60 


59.56 


60 


99.26 




(60.18) 


(0.61) 


(59.57) 




(99.27) 


130 


66.12 


0.54 


65.58 


66 


99.36 




(66.15) 


(0.56) 


(65.59) 




(99.38) 


132 


72.12 


0.51 


71.61 


72 


99.46 




(71.83) 


(0.46) 


(71.37) 




(99.13) 


134 


78.14 


0.48 


77.65 


78 


99.56 


136 


83.59 


0.47 


83.12 


84 


98.95 


138 


88.92 


0.44 


88.48 


90 


98.31 


140 


94.69 


0.65 


94.04 


96 


97.96 


142 


100.46 


0.43 


100.03 


102 


98.07 
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FIG. 1. HF(Sk3)+BCS energy of the isotopes 124_142 Xe as a function of the quadrupole deformation /3 for fixed pairing gap 
(dashed line) and for fixed pairing strength (solid line) treatments (with 1 MeV of separation between absolute minima), as 
well as deformation dependence of pairing gaps from the fixed strength calculation (dashed line for proton gap, dotted line for 
neutron gap). The scale in the vertical axis is 1 MeV between two ticks. Vertical lines indicate ground states from the fixed 
gap calculation. 
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FIG. 2. Discrete and gaussian-folded Gamow- Teller strength distributions B(GT _ ) for 124 ~ 142 Xe, from a HF(Sk3)+BCS(G) 
calculation. Two equilibrium deformations appear for 124_132 Xe, the solid line and the discrete spectrum corresponding to the 
prolate shape, and the dashed line corresponding to the oblate shape. The calculations for 134 ' 136 > 13S Xe correspond to spherical 
shapes, while those for 140 ' 142 Xe correspond to prolate shapes (see Table II). 
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FIG. 3. The same as in fig. 2 but for GT strength distributions. 
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FIG. 4. GT strength distributions in 128 Xe (upper panels) and 140 Xe (lower panels). Left panels show GT strengths, while 
right panels show GT + strengths. Calculations with prolate (spherical) shapes are plotted upward (downward), respectively. 
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FIG. 5. Gaussian-folded Gamow- Teller strength distributions B(GT ± ) in Pa/(4tt) units for both equilibrium shapes of 128 Xe, 
from a HF(Sk3)+BCS(G) calculation (solid line) and a HF(Sk3)+BCS(A) calculation (dashed line). 
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